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Abstract 

Derivative polynomials in two variables are defined by repeated differentiation of the 
tangent and secant functions. We establish the connections between the coefficients of these 
derivative polynomials and the numbers of interior and left peaks over the symmetric group. 
Properties of the generating functions for the numbers of interior and left peaks over the 
symmetric group, including recurrence relations, generating functions and real-rootedness, 
are studied. 
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1 Introduction 

Let S n denote the symmetric group of all permutations of [n], where [n] = {1,2, ... ,n}. A 
permutation ir = 7r(l)7r(2) • • • 7r(n) £ S n is alternating if 7r(l) > vr(2) < •••7r(n). In other 
words, ir(i) < ir(i + 1) if % is even and > w(i + 1) if i is odd. Similarly ir is reverse 
alternating if 7r(l) < 7r(2) > •••ir(n). Let E n denote the number of alternating permutations 
in S n . For instance, E4 = 5, corresponding to the permutations 2143, 3142, 3241, 4132 and 
4231. The number E n is called an Euler number because Euler considered the numbers E^n+i- 
The bijection tt h- > tt c on S n defined by vr c (i) = n + 1 — ir(i) shows that E n is also the number 
of reverse alternating permutations in S n . Alternating permutations have rich combinatorial 
structure and we refer the reader to the survey paper by Stanley [18] for recent progress on this 
subject. 

In 1879, Andre [1] observed that 



Ex 
E n — - = tan x + sec x, 
m 1 



71! 
n=0 



X 2 „X 3 „X 4 ^X h 
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Since the tangent is an odd function and the secant is an even function, we have 

00 x 2n+l 00 x 2n 

EE2n+i~, — =tanx and > Eo n -, — rr = secx. 
+ 2n+l ^ (2n) 

n=0 v 1 n=0 y ' 

For this reason the integers -E?2n+i are sometimes called the tangent numbers and the integers 
E2n are called the secant numbers. 

Let D x denote the differential operator -4-. In 1995, Hoffman [12] considered two sequences 
of derivative polynomials defined respectively by 

D x {taxi{x)) = P n (tan(x)) and D x (sec(x)) = sec(x)Q n (tan(x)) (1) 

for n > 0. From the chain rule it follows that the polynomials P n {u) satisfy Pq{u) = u and 
P n+ i(u) = (1 + u 2 )P' n {u), and similarly Qo(u) = 1 and Q n+ i(u) = (1 + u 2 )Q' n (u) + uQ n (u). 
In particular, the numbers P n (0) and Q n (0) are respectively the tangent and secant numbers. 
Recently, several combinatorial formulas for these derivative polynomials have been extensively 
investigated (see [81 [TTJ [131 E] f° r instance). For example, let 

tan fc (x) = Y T(n,k) — 
f n! 

n>k 

and 

^ x n 

sec(x) tan fc (x) = ^j5(n,A;) — -. 

n>k 

The numbers T(n, k) and S(n, k) are respectively called the tangent numbers of order k (see [U 
p. 428]) and the secant numbers of order k ((see [5j p. 305])). Clearly, the numbers T(n, 1) and 
S(n, 0) are the tangent and secant numbers, respectively. Cvijovic [H Theorem 1] obtained that 

n+l ^ 

P n {x) = T(n, 1) + ^T{n + 1, k)x k 
k=l 

and 

n 

Q n (x) = Y J S{n,k)x k . 

k=0 

The organization of this paper is as follows. In Section 2, derivative polynomials in two 
variables are defined by repeated differentiation of the tangent and secant functions. In Section 
3, we study the connections between interior peaks polynomials and left peak polynomials. In 
Section 4, we obtain several explicit formulas. In Section 5, we give both central and local limit 
theorems for the coefficients of certain polynomials. 

2 Triangular arrays 

Since D x (t&n(x)) = sec 2 (x) and D x (sec(x)) = tan(x) sec(x), it is convenient to set y = tan(x) 
and z = sec(x). Hence we have D x (y) = z 2 and D x (z) = yz. It is natural to consider the triangu- 
lar arrays (H / n,fc)„>i,o<fc<L(n-i)/2j , (W"i ifc )n>i,o<fc<Ln/2j and (R n ,k)n>i,o<k<n defined respectively 
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by 

L(n-1)/2J \n/2\ 

= E W^y^-'z 2 ^ 2 , D n x (z) = £ Wi tk y n - 2k z 2k+1 (2) 

fc=0 fc=0 

and 

n 

D^y + z) = Y,Rn, k y n - k z k+1 . (3) 

fc=0 

In the following discussion we will always assume that n > 1 and < A; < n. By the linearity 
of D x , we have 

W n fc-i if /c is odd, 



W k if A; is even. 



(4) 



Theorem 1. For n > 1 and < < n, we have 

Rn+i,k = {k + l)Rn,k + {n-k + 2)R n)k -2, (5) 
W n , k = (2k + 2)W n _i ik + (n - 2A;)VF n _ lifc _ 1 (6) 

and 

Wi >k = (2k + l)W l n _ 1>k + (n-2k + l)W l n _ 1)k _ v (7) 

Proof. Note that 

n n 

+ z) = D x (D n x (y + z)) = + l)R n , k y n - k+l z k+1 + J> " k)R n ^ k y n ~ k ~ 1 ^ 



y k+3 

10 — nj±t r ^ k y 

k=0 k=0 



Thus we obtain ([5]). Similarly, we get ([6]) and ([7]). □ 

The numbers W nk and fe arise often in combinatorics and other branches of mathematics 
(see H3 A008303, A008971]). In US US], Petersen studied the peak statistic over S n . Let 
ir = 7r(l)7r(2) • • • 7r(n) G S n . An interior peak in 7r is an index i G {2, 3, . . . , n — 1} such that 
7r(i — 1) < 7r(i) > ir(i + 1). Let pk (ir) denote the number of interior peaks in ir. An left peak in 7r 
is an index i £ [n — 1] such that 7r(i — 1) < ir(i) > ir(i + 1), where we take 7r(0) = 0. Let pk' (vr) 
denote the number of left peaks in ir. For example, the permutation ir = 21435 has pk (ir) = 1 
and pk' (vr) = 2. 

Let 

W n (x) = E z pkW W l n (x) = E ^ pk ' W - 

Note that degW n (x) = - 1)/2J and degW^(x) = [n/2\. Using © and ©, it is easy to 
verify that W nj jt is the number of permutations in S n with k interior peaks, and W l n k is the 
number of permutations in S n with k left peaks. Hence 

W n (x) = Y,W n , k x k and ^J^^. 

fe>0 fe>0 
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Therefore, by Q, the polynomials W n (x) satisfy 

W n+1 (x) = {nx-x + 2)W n (x) + 2x(l - x)D x W n (x), 

with initial values W\(x) = 1, W^(x) = 2 and Ws(x) = 4 + 2x. By ©, the polynomials W^(x) 
satisfy 

W^s) = (rax + l)W l n {x) + 2x(l - x)D x W l n {x), 

with initial values W-j^x) = 1, W^ix) = 1 + x and W^(x) = 1 + 5x. 

We now present a connection between the polynomials W n (x) and W^(x). 

Theorem 2. For n > 1, toe /iave 

W n (u) = (U - l) 2 * 1 U- 1 P n ((u - 

and 

W^(«) = («-l)5Q fl ((u-l)-3). 

Proof. Set y = tan(x) and z = sec(x). Note that z 2 = y 2 + 1. Comparing (pQ) with ([2]), we have 

L(n-1)/2J 

Pn(y) = £ W^ fc y n - 2fc - 1 (l + y 2 ) fc+1 = (y^ 1 + y n+1 )T^(l + y~ 2 ), 

fc=0 

and 

[n/2\ 



Qn(y) = £ < fc y n " 2fc (i + y 2 ) fc = y n ^(i + y" 2 )- 

k=0 

Let u = 1 + y -2 . Then the desired results follows immediately. □ 

Let R n {x) = Yl=Q R n,kX k . By ©, we have 

R n+1 (x) = (l + nx 2 )R n (x)+x(l-x 2 )R' n (x) for n > 1, (8) 

with initial values Ri(x) = 1 + x, R2(x) = 1 + 2x + x 2 and i?3(x) = 1 + 4x + 5x 2 + 2x 3 . Set 
i?o(x) = 1. For < n < 6, the coefficients of R n {x) can be arranged as follows with R n< k in row 
n and column fe: 

1 

1 1 

1 2 1 

14 5 2 

1 8 18 16 5 

1 16 58 88 61 16 

1 32 179 416 479 272 61 

Using (dD, we get 

R n (x) = xW n {x 2 ) + W l n {x 2 ). (9) 
Therefore, R n $ = 1, R nil = 2™" 1 , Y2=o R n,k = 2n! and i? n , n = iV-i,n-2 = -En for n > 2. 
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3 Identities 



Recently, much attention has been paid to the properties of the polynomials W n (x) and W^(x). 
For a permutation tt E S n , we define a descent to be a position i such that > n(i + 1). 
Denote by des (n) the number of descents of ir. The generating function for descents is 

A n (x)= * desW > 

7TS>Sn 

which is well known as the Eulerian polynomial. The exponential generating function for A n (x) 
is 

At \ , , ,2" (1 - x)e Z ^- X ^ , s 

A(x,z) = 1 + £ A(*)^ = 1 _ J eZil - x) • (10) 

n>l 

By the theory of enriched F '-partitions, Stembridge [HI Remark 4.8] showed that 

In [151 Observation 3.1.2], Petersen observed that 
4x 



) = (T^T ((! " *)" + £ ( ■) (1 - *)" ^'■'■M-r^j • (12) 



n V(i + ^)V ( 

In order to simplify the identities (jlip and (|12p . we will present another connection between 
W n (x) and W^(x). Let C n denote the set of signed permutations of ±[n] such that oj(— i) = —uj(i) 
for all i, where ±[n] = {±1, ±2, . . . , ±n}. Let 



^.ades (w) 

where 

des (cj) = \{i e [0, n - 1] : > w(i + 1)}| 

and 

ades (ijj) = \{i G [0, n] : > uo(i + 1)}| 
with w(0) = cj(n + 1) = 0. For example, if 

oj = (-2,-4,6,-8,1,3,7,5), 

then des (ui) = |{0, 1, 3, 7}| =4 and ades (cj) = |{0, 1, 3, 7, 8}| = 5. Recently, Dilks, Petersen and 
Stembridge [1U| studied the combinatorial expansions for the affine Eulerian polynomials C n (x). 
From [TOj Cor. 5.6 and Cor. 5.7], we get 



2xW n 

and 

{l + x)W, t 



Ax \ C n (x) 



(l + x)V (1 + x)™" 1 
4x \ C n (x) 



(l + x) 2 J (1 + x)"- 1 ' 
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For n > 1, set 

T n {x) = C n {x 2 ) + -C n {x 2 ). 

x 

From their generating functions (see [TOj Sec. 6.2]), we get 

C(X, Z) = l + J2 d n(x)^ = - _ X ~/ Z[x _ x) , 
n>l 

C(x, *) = ! + (! + x)z + ^ C n {x)- } = \ _ J . 

n>2 

Therefore, 

T(x, z) = l + J2 T n {x) Z - = C(x\ z) + -(C{x\ z)-l) = 1 — * (13) 

^-^ n\ x 1 — xe z " x > 

n>l 

Comparing (fTU|) with (fT3"j) . we get 

x + T(x,z) = (l + x)A(x,z(l + x)). (14) 
Using (|14p . we immediately obtain the following result. 
Theorem 3. For n > 1, we have 

T n (x) = (l + x) n+1 A n (x). (15) 

Motivated by the connections between W n (x) and W^(x), in the following sections we will 
study the properties of the polynomials R n (x) . 

4 Explicit formulas 

The goal of this section is to find an explicit formula for the polynomials R n {x). Let 

z n 

P(x,z) = TR n (x) — . 

^-^ n\ 

n>0 

Multiplying both sides of (JH]) by x n jn\ and summing over all values of n, we get that P(x, z) 
satisfies the following partial differential equation: 

x{x 2 _ 1} ^<M + (1 _ x 2 z) ^A = z) + x . (16) 

OX OZ 

It is well known pU A008303, A008971] that 

n 

W(x,z) = Y^W n (x) — 

sinh(z\/l — x ) 



n>l 



\Jl — xcosh(2\/l ~~ x ) ~ sinh(z\/r" 



and 



n! 

n>l 



\/T — xcosh(z\/l ~~ x ) ~ sinh(z\/l — x) 

Using (|9|), we get 



P(ac, z) = 1 + V [*(i 2 ) + w'fac 2 )]^ = xW(x 2 , z) + W'(x 2 » z). 

n>l 

Thus 

xsinh(z\/l - £ 2 ) + Vl - x 2 

P(x, z) = . = ^=^= , (17) 

V 1 — x 2 cosh (zv 1 — x 2 ) — sinh(^v 1 — x 2 ) 

and this formula gives a solution to the partial differential equation (|16jl . It should be noted 

that the generating function of the polynomials R n {x) has also been studied by Charalambos [6l 

p. 542]. Recall that the hyperbolic cosine and the inverse of the hyperbolic cosine are defined 

respectively by 



e + e i 

coshz = and arccosh (x) = ln(x + y x 2 — 1) for x > 1. 

Let 

t n 

R(x,t) = ^ R n+1 (x) — . 

^— ' n! 

n>0 

Clearly, 

^ =*.,«). 

By i/ie method of characteristics (see [20] for instance), Charalambos [BJ p. 542] obtained a 
two-variable generating function 

""■" ^l"f-D ' (18) 



with z = —t\J 1 — x 2 + arccosh (1/x). Comparing (JT7J) with f)18[) . the latter is simpler. To obtain 
an explicit formula for R n (x), we will use (|18H . Let {xj}j>i be a sequence of variables. The 
partial Bell polynomials B n ^ =: B n ^(x\,X2, ■ ■ ■ ,x n -k+i) are defined by the generating function 

n>k \i>l / 

with i?o,o = 1 and B n> o = for n > (see [3 p. 133]). 

Proposition 4. Lei P ni £ be the partial Bell polynomials. When Xi = (l — x 2 ) 1 ^ 2 "' /or eac/i 
i>l, we have 

n 

R n+l (x) = ^(-lf-^Kl + x) fc+1 P n ,*- (20) 
fc=i 
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Proof. Let z = —tyl — x 2 + arccosh (1/x). Note that 



_|_ g— z g—t^l—x' 2 _|_ gt\/l— x' 2 g— ty/T— x 2 gi\/l— x' 2 

x( - 1) = - + \A -x 2 - X 

where yi = (l — x 2 ) L( l+1 )/ 2 J _ Thus 

i?(x,t) 



i-^ + Ei>i(-i)^| 
1 

l + x 



i + (i + ^)Ei>i(-i)^S' 



where xi = (l — x 2 ) ^ ^ . Then the desired result follows immediately from the formula for 
the geometric series and (I19p . □ 



Here we provide an example for illustration of Proposition SJ 
Example 5. Consider the case n = 4, we have 

B^i = £4, i?4 5 2 = 4xiX3 + 3X 2 ,, -^4,3 = 6x 2 X2, -64,4 = X 4 

fsee |7| p. 507/ /or instance). When Xi = (l — x 2 ) * /or eac/i i > 1, we <?e£ 

B 4 ,i = 1 - x 2 , 5 4 ,2 = 7 - 4x 2 , £4,3 = 6, £4,4 = 1. 

4 

^(-l) 4 ~ fc /c!(l + x) k+1 B i>k = 1 + 16x + 58x 2 + 88x 3 + 61x 4 + 16x 5 . 
fc=i 

Recall that 

S n ,fc(l, 1,1,...) = S(n,k), 

where S(n,k) is a Stirling number of the second kind (see [71 p. 135]). Hence when x = 0, the 
formula (l20l) reduces to 



uQ-ir^!%fc). 

k=0 

Moreover, when x = 1, the formula ([20]) reduces to 

/? 

(n + 1)! = ^(-l) n - fc A:!2 fc J B nifc (l, 1, 0, 0, ... , 0). (21) 
fc=i 

From (|20p . we observe that R n {x) is divisible by (1 + x) 2 for n > 2. It is well known that 
the Eulerian polynomial A n (x) has only real nonpositive simple zeros (see [31 Theorem 1.33]). 
In the next Section, we will show that the polynomial R n (x) also has only real zeros. 
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5 Central and local limit theorems 



Let sgn denote the sign function denned on IR, i.e., 



sgnx 



+1 if x > 0, 
if x = 0, 
-1 ifx<0. 



Let RZ denote the set of real polynomials with only real zeros. Furthermore, denote by RZ(/) 
the set of such polynomials all whose zeros are in the interval /. Suppose that f,F£ RZ. Let 
{r{\ and {sj} be all zeros of / and F in nonincreasing order respectively. We say that / separates 
F, denoted by / < F, if deg / < deg F < deg / + 1 and 

si > T\ > s 2 > r 2 > s 3 > r 3 > • • • . 

It is well known that if / G RZ, then /' G RZ and f < f. 

Theorem 6. For n > 1, we have R n {x) £ RZ[— 1,0) and R n {x) ^ R n +i(x). More precisely, 
R n (x) has |~n/2] — 1 simple zeros, and the zero x = — 1 with the multiplicity [n/2\ + 1. 

Proof. From flSJ), we have deg R n+ i(x) = degR n (x) + l and R n (0) = 1. Note that Ri{x) = 1 + x, 
R 2 (x) = (1 + x) 2 , R 3 (x) = (1 + x) 2 (l + 2x) and i? 4 (x) = (1 + x) 3 (l + 5x). So it suffices to 
consider the case n > 4. We proceed by induction on n. Let m = [n/2\ + 1. Suppose that 



Rn(x)=H(x-r l )(x + l) m , 



i=l 

where > r\ > r% > ■ ■ ■ > rf. > — 1. 

Let g(x) = Ui=i( x ~ n) and G(x) = R n+1 (x)/(x + l) m . Then degG(x) = deg^(x) + 1. 
By (JHJ) , we have 

G{x) = (1 + nx 2 )g(x) + x(l - x 2 ) g ^ XJ + toc(1 - x)5(a;). 

X — Tn 

Note that sgnG(rj) = (-1)' for 1 < i < k and G(-l) = (1 + n - 2m)g{-\). Hence G{x) has 
precisely one zero in each of k intervals (r^, rk-i), ■ ■ ■ , (?"2, ri), (ri, 0). Note that 



1 + n - 2m = n - 2[n/2j - 1 = < 



if n is odd; 
— 1 if n is even. 



Thus —1 is a simple zero of G(x) if n is odd. If n is even, then sgn G(—l) = (— and G(x) has 
precisely one zero in the interval (— 1, r^). Hence i? n+ i(:r) G RZ[— 1,0) and R n (x) ^ i? n+ i(x). 
This completes the proof. □ 
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For n > 1, set 

R n (x) = (l + x)^ +1 G n (x). 
Then the polynomial G n (x) has only positive integer coefficients. 



(22) 



Open Problem 7. Is the similarity between \15)) and A22\) just a coincidence? Can either 
equation be given a combinatorial interpretation? 

Let {a(n,k)}o<k<n be a sequence of positive real numbers. It has no internal zeros if there 
exist no indices i < j < k with a(n,i)a(n,k) ^ but a(n,j) = 0. Let A n = X]fc=o a ( n '^)- 
We say the sequence {a(n,k)} satisfies a central limit theorem with mean /i n and variance o\ 
provided 



lim sup 



£ 

fc=0 



a(n, fe) 

-An. 



1 



% 2 

-oo 



0. 



The sequence satisfies a local limit theorem on B 6 K if 

<7 n a(n, /i n + xcr n ) 1 



lim sup 



Recall the following Bender's theorem. 



2ir 



0. 



Theorem 8. fM/ Let {-P n }n>i be a sequence of polynomials with only real zeros. The sequence 
of the coefficients of P n satisfies a central limit theorem with 

2 



li n = -^IJ. an d at = J 1 )/, + 



p n (l) ' p n (l) \p n (l) 



provided that lim az = +oo. If the sequence of coefficients of each P n (x) has no internal zeros, 
n— >oo 

then the sequence of coefficients satisfies a local limit theorem. 

Combining Theorem [6] and Theorem [81 we obtain the following result. 

Theorem 9. The sequence {R n ,k}o<k<n satisfies both a central limit theorem and a local local 
limit theorem with fi n = (2n — l)/3 and a\ = (8n + 8)/45, where n > 4. 

Proof. By differentiating 0, we obtain the recurrence 

x n+ i = (4n)n! + (n - l)x n 

for x n = R' n (l), and this has the solution x n = (An — 2)n!/3 for n > 2. By Theorem [81 we have 
/i n = (2n — l)/3. Another differentiation leads to the recurrence 

?/n+i = ^n!(4n 2 - 5n + 3) + (n - 3)y n 

for y n = -R'^(l). Set y n = (an 2 + bn + c)n\ and solve for a, 6, c to get y n = n\(A0n 2 — 84n + 56)/45 
for n > A. Hence a\ = (8n + 8)/45, where n > A. Thus lim a 2 = +oo as desired. □ 
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Let P{x) = Y^i=o a i xt be a polynomial. Let m be an index such that a m = maxo<j< n aj. 
Darroch [9] showed that if P(x) S RZ(— oo,0], then 



So the following result is immediate. 

Corollary 10. Let i be an index such that i? n/ j = maxo<fc< n Rnk- If ~ l)/3 is an integer, 
then i = (2n- l)/3; Otherwise, i = [(2n - l)J/3 or i = \{2n - l)]/3. 
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